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ABSTRACT 

We propose a mechanism of long-range coherent coupling between nuclear spins to be 
used as qubits in solid-state semiconductor-heterojunction quantum information processing 
devices. The coupling is via localized donor electrons which in turn interact with the two- 
dimensional electron gas. An effective two-spin interaction Hamiltonian is derived and 
the coupling strength is evaluated. We also discuss mechanisms of qubit decoherence and 
consider possibilities for gate control of the interaction between neighboring qubits. The 
resulting quantum computing scheme retains all the gate-control and measurement aspects 
of earlier approaches, but allows qubit spacing at distances of order 100 nm, attainable with 
the present-day semiconductor device technologies. 



PACS: 73.20.Dx, 71.70.Ej, 03.67.Lx, 76.60.-k 



Recent technological advances in electronics related to spin polarization [1,2] have 
boosted experimental and theoretical interest in quantum information science in condensed 
matter systems, specifically, in semiconductor heterostructures at low temperatures and 
in high magnetic fields. The solid state implementations of quantum information devices 
seem to be among the most promising ones, due to possible scalability of the elementary 
logic gates into more complicated integrated circuits. Several designs for solid state and 
related spin-based quantum information processors have been suggested [3-8] . Preliminary 
experiments, involving several quantum bits (qubits), have been carried out or are being 
contemplated [9-10]. 

Our work stems from the proposals that utilize nuclear or electronic spins as qubits 
for information processing [3-7]. These are natural choices for qubits because at low tem- 
peratures spin states in semiconductors have relatively long decoherence times, sometimes 
milliseconds or even longer for electronic spins, and seconds for nuclear spins [11-14]. We 
propose a new mechanism for coupling between two nuclear-spin qubits, combining aspects 
of two models of quantum information processors, one based on nuclear spins in quantum 
Hall effect systems [4], and another utilizing the nuclear spins of phosphorus donors in a 
silicon heterostructure [5] . 

An appealing aspect of Kane's model [5] is a possibly experimentally feasible scheme 
for reading out the state of the quantum register, i.e., measurement of a nuclear spin, 
achieved by transferring the nuclear spin polarization to the electronic state, while the later 
is measured with the use of a single electron transistor. The model proposed in [4] has a 
different advantage: unlike [5] , the interaction between the nuclear spins is mediated by the 
two-dimensional (2D) electron gas, and thus is longer ranged due to the highly correlated 
state of the 2D electron gas in the quantum Hall regime. This opens up possibilities for 
experimental realization of such quantum information processors, because large separation 
between spin qubits means greater lithographic dimensions in manufacturing the device. 
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The price paid is that the coupling is weak, and therefore the time scales of the "gate 
function" can be as large as 1 s. 

In this work we combine the two proposals, thus retaining the measurement and 
control scheme proposed in [5,7,9] and at the same time allowing larger separations, of 
order 100 nm, between interacting qubits. The resulting system is thus realizable with the 
present-day semiconductor technologies. We propose a model where sparsely positioned 
phosphorus donors are imbedded in a 2D electron gas in the quantum-Hall regime. The 
localized donor electrons interact via the delocalized 2D electrons and thus indirectly me- 
diate nuclear-spin interactions. In 3D, spin coupling mechanisms via conduction electrons 
have been well studied [15]. Here we estimate the range of this induced nuclear-spin in- 
teraction for the 2D case and find it to be of order 100 nm. This is large compared to 
atomic dimensions, donor-electron bound state radii, and even the electronic magnetic 
length which is typically of order 10 nm. We find that this interaction is also stronger, 
thus corresponding to faster "gate function" times, than in [4]. 

We assume that the coupling between the electronic and nuclear donor spins is given 
by the Fermi contact interaction, H e - n = Aa n ■ a e . Here A = (87r/3) / u_e(7 n/ u n |\l/o(0)| 2 , 
where fx n and g n are the nuclear magneton and nuclear g- factor, respectively, |\I/o(0)| 2 
is the donor electron probability density at the nucleus, \ib denotes the Bohr magneton, 
and cr's are Pauli matrices. Coupling of the delocalized electrons to the nuclear spin is 
considerably weaker than that of the localized donor electron. Therefore, we assume that 
the nuclear spin interacts with conduction electrons indirectly via the donor electron. 

As a prototype system, we consider 31 P donors positioned in Si, so all the spins in- 
volved are 1/2. The donor electronic and nuclear spins form a four-level system. The spec- 
trum of this two-spin system can be obtained to O(A) with H e - n treated as perturbation. 
The energy levels are E = -( 7n + A)/2 + A, E x = ( 7n -A)/2-A, E 2 = (- 7n + A)/2-A, 
and E3 = ( 7n + A)/2 + A, where 7n = g n fi n H is the nuclear spin splitting. Here 
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H is the magnetic field, and the expression for A, the electronic Zeeman gap, will be 
given shortly. The eigenstates associated with these energy levels are |0) = | | e | n ), 



1) = I Un> + (24/ A) | TeU), |2> = | UU) - (2A/A)\ |etn), and |3> = I Tetn), where 



I leln) represents the electronic and nuclear spin down state, etc. Here we propose to 
consider the states |0) and |1) as qubit states of a quantum computer. By altering the 
hyperfine coupling constant A by distorting the spatial state ^ of the donor electron with 
an electrostatic gate [5,7], one can selectively control the state of an individual qubit by 
means of the NMR technique. 

In order to calculate the interaction Hamiltonian between two qubits, we first consider 
the coupling between the donor electron and conduction electrons. The ground state of the 
donor electron is bound (localized) and will typically lie in the energy gap, several meV be- 
low the conduction band edge. For temperatures of order mK, electronic transitions from 
this localized state to the conduction band are highly improbable. The dominant interac- 
tion between the localized electron and conduction electrons is their Coulomb interaction. 
We are interested only in the exchange part of this interaction, i.e., the spin-dependent 
part. The spin-independent part causes screening, but it is weak in 2D [16] and, especially 
in the presence of the magnetic field, cannot ionize the donor. 

In large magnetic field, the delocalized 2D electrons occupy highly degenerate Landau 
energy levels [16]. It is convenient to introduce electron bound state creation and anni- 
hilation operators and 6 ns , where n represents the donor spatial state, and s is the 
spin ^-component, j or {. Let a mk s , a m k xS denote the creation and annihilation opera- 
tors for the delocalized 2D electrons, where m labels the Landau level, while Tik x is the 
x-momentum (we use the asymmetric gauge). Then the exchange coupling between the 
bound and delocalized electrons can be written as 



where the sum is over all the indices. Here we have neglected the spin-orbit interaction. 
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(1) 



In what follows, we will retain only the lowest donor electron spatial state, i.e., account 
only for the transitions between the two Zeeman levels of the ground state. 

The 2D electrons are assumed to be in a nondissipative quantum Hall state with filling 
factor v = 1, i.e., the lower Zeeman sublevel of the Landau ground state is completely filled 
[4] . This choice ensures reduced decoherence and relaxation effects [14] , owing to the energy 
gap in the spectrum of the lowest-energy spin-wave excitations which are well studied 
[17,18]; their spectrum is given by £ k = A + E c [l - I (£ 2 k 2 /4) exp (-£ 2 k 2 /4)] , where I 
is the modified Bessel function. Here A = g^sH is the Zeeman gap, E c = (n/2) 1 ^ 2 {e 2 /t£) 
is the characteristic Coulomb energy, and g is the effective g-factor in the potential well 
that holds the 2D electron gas, while e is the dielectric constant of the material, and 
£ = (hc/eH) 1 / 2 is the magnetic length. Extension to larger integer filling factors is possible 
[14,17,18]. One can also introduce [18] normalized creation and annihilation operators for 
the spin waves, quadratic in electronic operators, 

St-^Y^V^W a k (2) 

Here L XyV are the transverse dimensions, taken to infinity in the final calculation. The 
summation over p is taken in such a way [18] that the wave number subscripts are quantized 
in multiples of 2n/L x . The spectrum of these spin- waves has been experimentally verified 
in GaAs heterostructures [19]. 

We will include only these lowest excitations in the sum (1); our goal is to rewrite 
(1) in terms of the spin- wave operators (2). The exchange coupling is thus truncated to 

G ^Cm',k x ,k' x = G kx ,k'Jn,oS n ',oSm,oSm',o, where 

G KK = J rf 3 R 1 rf 3 R 2 vl/S(R 1 )vl/o(R2)^(Ri - R 2 )^,^(R 2 )$o,^(Ri), (3) 

U(Ri — R2) = e 2 /e|Ri — R 2 | is the Coulomb interaction and \I/o(R) is the donor electron 
ground state. The states of the conduction electrons confined in the 2D well are $0,/^ (R) = 
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(f>o : k' x ( r )x( z )i where 0o,fc^( r ) are the standard 2D Landau states [16]; x( z ) describes the 
confinement of the conduction electron wavefunction in the z direction and depends on 
the nature of the confinement potential. Here and in the following R = (r, z) , with R and 
r = (x, y) being 3D and 2D coordinates, respectively, while z is the direction perpendicular 
to the heterostructure, in which the applied magnetic field is pointing. 

With the use of the expressions for Landau ground state wavefunctions, 4>o,k' x (r) = 
£~ 1 (nL x )~ 1 / 2 e' lkxX exp[— (y — £ 2 k x ) 2 /2£ 2 ], and (2), after a lengthy calculation, we get 

H ex =\Y.\ W *\ teXie |Sk + W£\ jeXTe |^] (4) 
k 

where | te)(le I = &}&| in the appropriate subspace, and 
Wk = £(2nL x L y )V 2 J rf3R i rf3R 2^o(Ri)^o(R 2 )^(Ri - R 2 )x*(* 2 )x(*i)Ck(ri, r 2 ) (5) 



C k (ri, r 2 ) = exp [(xi ~ x 2 ) 2 + (yi - y 2 f - 2i(xi - x 2 )(yi + 2/2)] | 



x 



exp 



£ 2 k k 

-{kl + k 2 ) - -^-(iyi + iy 2 - x\ + x 2 ) - -§-{ix\ + ix 2 +yi- y 2 ) 



(6) 



Note that since all the position vectors R, r are measured from the origin at the donor 
atom, the quantity Wk depends also on the donor coordinates. To the leading order, (4) 
gives the interaction of the donor electron spin with excitations of the 2D electron gas in 
the v = 1 integer quantum-Hall state. 

One can rewrite the interaction (4)-(6), with (4) multiplied by the unit operator in 
the nuclear-spin Hilbert space, in terms of the eigenstates of the electron-nucleus system. 
With the use of the expressions derived earlier for these eigenstates in terms of direct 
products of electronic and nuclear spin states, we obtain 

H ex = 2 E w * (x |1)(01 + 12)(01 + 13)(11 - x |3)(2| J 5k + h - c - (7) 



Now one can calculate an effective Hamiltonian for the interaction of two qubits. Since 
the electronic Zeeman gap is much larger than the nuclear one, we can truncate the Hilbert 
space of the combined electron-nucleus spins to the two lowest lying states. Thus we retain 
only the |0)(1| and conjugate transitions in the exchange interaction (7). 

An effective interaction between two qubits can be obtained within the standard frame- 
work of second order perturbation theory by tracing out the states of the spin waves; see 
[15,20,21] for similar calculations. The result can be written as 



The subscripts 1 and 2 in (8)-(9) label the two donor qubits, while Wk,i and W^z are the 
coupling constants of each donor electron spin to spin waves, given by (5), and £k is the 
spin- wave energy. 

The nuclear-spin energy gap is much smaller than the electronic spin-wave excitation 
energies. Therefore, we can ignore E\ — E in the denominator in (9). Furthermore, due to 
the large value of the spin- wave spectral gap at k = 0, £q = A, we do not have the "small 
denominator" problem encountered in other calculations of this sort, e.g., [20]. Physically 
this means that the spin excitations in the 2D electron gas mediating the effective qubit- 
qubit interaction are virtual and so this interaction does not cause appreciable relaxation 
or decoherence on the "gate function" timescale %/J. 

It is important to note that one can construct a universal CNOT logic gate from 
the controlled dynamics governed by Hamiltonians of the form of H\ i and single qubit 
rotations [6]. The coupling strength J between the qubits can be externally controlled 
by the electrostatic gates built above the 2D inversion layer. By applying gate voltages, 



Hi >2 = J|0il 2 )(li0 2 | + J*|li0 2 >(0il 2 



(8) 



Here the coupling constant between the two qubits is 
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one can locally vary the density of the 2D electrons thus changing coupling between the 
delocalized and donor electrons. This results in control over the effective coupling constant 
J in (9). The precise effect of gates on interactions between the qubits as well as on 
decoherence of their states, should be further studied in order to establish the feasibility 
of the quantum-computing approach proposed here. Most other semiconductor solid-state 
quantum-computing approaches [3-7] utilize gates. 

Let us explicitly calculate the coupling constant J in (8)-(9). Because the spatial 
ground state of the donor is localized on a scale smaller than the magnetic length £, 
the overlap integrand in (5) is vanishingly small for |ri — r 2 | > £. At the same time, 
for |k| > l/£, the value of Ck decreases exponentially. Thus Ck can be simplified by 
neglecting x\ — x 2 and y\ — y<i terms in (6). Moreover, for two donors at separation 
larger than £, we can put (ri + r 2 )/2 ~ r J7 with r,- being the location of either one of 

them. Then (5) can be approximated by Wkj = ^(L x L y ) -1 / 2 exp ^— ^-j «k • r^j, with 

Z= (l/27r£ 2 ) 1/2 /rf 3 Ri^ 3 R-2^(Ri)^o(R2)^(Ri -H 2 )x*(z2)x(zi)- 

Finally, the coupling constant J of the effective interaction (8) can be obtained by 
transforming the summation in (9) to integration in the limit L XyV — > oo, 

fe(;) 1/2exp (-s)' < r >^ < 10 > 

where d = (E C /2A) 1/Z £. A similar dependence of the coupling on the donor separation 
r was obtained in a study of nuclear polarization diffusion in the quantum Hall regime 
[21]. Interaction (7)-(8) between the spins has finite range d, which, however, is very 
large compared to the effective Bohr radius of the donor ground state. Thus the indirect 
exchange at large distances dominates the direct exchange interaction resulting from the 
overlap of the two atomic wavefunctions. For magnetic field H = 6T and e = 12, we get 
d — 65 nm, which is indeed much greater than the characteristic Bohr radius for a donor 
electron in silicon. 
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In order to estimate J, we have to evaluate the overlap integral Z. For an order-of- 
magnitude estimate we will assume that x( z ) * s constant inside the well and zero outside. 
Then Z ~ (2tt)- 1 / 2 (^)- 1 / rf 3 R 1 d 3 R 2 *o( R i)*o(R2)^(Ri - R 2 ), where S is the width of 
the well. We put 5 ~ 4nm. For \l/o(R), the donor ground state, we choose a spherically 
symmetric Hydrogen-like ground state with the effective Bohr radius as — 2nm. This 
is, of course, not the case in a realistic situation [22]. The ground state of the donor will 
be influenced by the band structure, by the magnetic field and by the confining 2D well 
potential, while the states of the conducting electrons will be distorted by the impurity 
potentials. We are not aware of a thorough study of these effects for our system. For the 
purposes of an order of magnitude estimate, however, a spherical state should be sufficient. 

Evaluating the integral for the Coulomb potential U, we obtain Z ~ (5a 2 B /165)E c . 
Assuming that the two donors are separated by the distance r = 100 nm and using the 
value 2A/h = 58 MHz from [4], we obtain the estimate J/h ~ 10 2 s _1 . 

The clock speed of the information processor just described appears to be a fraction 
of kHz and should be compared with the time scales for relaxation and decoherence. The 
leading mechanism for these at low temperatures is through interaction with impurities. 
It has been found theoretically [12,23] and confirmed experimentally [2] that nuclear spin 
relaxation in the quantum Hall regime is slow and strongly dependent on the impurity 
potentials; typically, the relaxation time Ti is of order 10 3 s. In our case the interaction of 
a qubit with the 2D gas is stronger, and as a result, the relaxation is expected to be faster. 
An estimate from formulas in [12,23] gives T\ ~ 1 s. There is, however, another important 
issue — decoherence, on time scales T 2 . Recently, this quantity has been calculated in the 
same framework, that is, when the interaction of the conduction electrons with impurities 
is taken into account [14]. The results of [14] can be adjusted for the present case and 
yield the estimate T 2 ~ 10 _1 s. 

The existing quantum error correction protocols require the quality factor, equal the 
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ratio of the gate-function clock time to decoherence time, not to exceed 10~ 5 [24]. Our 
estimates indicate that this is not the case for the present system. Actually, no quantum 
computing proposal to date, scalable by other criteria, satisfies this 10~ 5 quality- factor cri- 
terion. The values range from 10 _1 to 10~ 3 . The resolution could come from development 
of better error-correction algorithms or from improving the physical system to obtain a 
better quality factor. In our estimate of the decoherence time scale, we used parameters 
typical of a standard, "dirty" heterostructure with large spatial fluctuations of the impu- 
rity potential. These heterostructures have been suitable for standard experiments because 
they provide wider quantum-Hall plateaus. Much cleaner, ultra-high mobility structures 
can be obtained by placing the ionized impurity layer at a larger distance from the 2D gas 
or by injecting conduction electrons into the heterostructure by other means. 

Thus, our present quantum-computing proposal offers a clear direction for exploring a 
physical, rather than algorithmic, resolution to the quality factor problem. This possibil- 
ity should be further examined both experimentally and theoretically. Our new quantum 
computing paradigm suggests several interesting avenues for research. The effect of gates 
on the switching of qubit interactions and on decoherence requires further investigation. 
The first experimental realizations will probably involve only few qubits. The interac- 
tions of these may be significantly affected by the geometry, specifically, the edges, of the 
heterostructure. 
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